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Abstract 



We use Moser's normal forms to study chaotic motion in two-degree hamil- 
tonian systems near a saddle point. Besides being convergent, they provide a 
suitable description of the cylindrical topology of the chaotic flow in that vicin- 
ity. Both aspects combined allowed a precise computation of the homoclinic 
interaction of stable and unstable manifolds in the full phase space, rather 
than just the Poincare section. The formalism was applied to the Henon- 
Heiles hamiltonian, producing strong evidence that the region of convergence 
of these normal forms extends over that originally established by Moser. 
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I. INTRODUCTION 



Normal Forms (NF) are among the successful methods for either analytic or numeric 
studies of dynamical systems; by performing a suitable coordinate transformation, we even- 
tually obtain a more simple or dynamically "transparent" version of the original system. This 
approach can be formulated either for generic systems including two-dimensional maps [I[ 
or for hamiltonian systems || and conservative maps 0]. So, even when a given NF does not 
converge, due to small denominators or exact ressonances, it is of importance for numerical 
purposes. This is just what occurs when the NF is obtained around a stable point or orbit. 
Then, in spite of the well known divergence in this case, a truncation allows us to follow 
the motion for a long time with great precision. In fact, the nonconvergent case is the most 
considered in the literature [f|-[(| . 

The present work concerns the NF around a unstable point or orbit of a conservative 
system. For that, Moser demonstrated their convergence for both maps []7| and hamiltonian 
systems ||. Although convergent, this case did not receive much attention until recently, 
presumably because the particle remains a very short time in that region. Nevertheless, we 
will see that the Moser normal forms (MNF) are both convergent and a powerfull tool for 
the search for the basic structures of the chaotic motion rather than just following a specific 
orbit for a long time. 

The usefulness of the MNF for the study of conservative chaotic maps is already known 
in the literature. It allowed precise analytical computations of homoclinic points and the 
periodic points with long period, which accumulate in the homoclinic ones ||,[10|- Additional 



good numerical results were obtained even if small dissipative perturbations were added [11 



Area-preserving maps are usually only simplified reductions (appropriate Poincare sec- 
tions) of autonomous hamiltonian systems of two degrees of freedom. In particular, the very 



complex two-dimensional homoclinic tangle |L2| is already a reduction of the much more 
involved chaotic motion in the full phase space. However, extending the use of the MNF to 
the hamiltonian case allows us to study directly the proper four-dimensional chaotic flow. 
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This has previously been attempted in the literature [13| , but without taking full advantage 
of the method that we shall develop here. Another use of MNFs is found in |TJj] and takes 
advantage of its convergence to ascertain stability transitions of families of periodic orbits 
near hamiltonians' saddle points. 

In section |TI] we will develop the MNF approach for the case of a generic autonomous 
hamiltonian of two degrees of freedom around a saddle point, encompassed by Moser's 
convergence proof. In that vicinity, the flow's topology is cylindrical rather than toroidal, 
in the case of a chaotic regime |15|]. Firstly, we construct the relations linking the original 
system to the corresponding more "transparent" normalized system. This is done through 
a near identity polynomial coordinate transformation. We will see that, besides convergent, 
that transformation also reveals in a natural way, the cylindrical character of the topology. 
Both features turn the MNF into a powerfull tool. So, it was possible to compute, precisely 
for the first time, the continuous structures in the full phase space, underlying the homoclinic 
tangle in a Poincare section: the (un)stable manifolds which originate at the saddle point 
and at each neighbouring unstable periodic orbit, the homoclinic orbits associated with the 
latter and the periodic orbits with long period which accumulate on the homoclinic orbits. 
In section |TJ we obtain the recurrence relations for the coefficients involved in the theory. 

In sections fV| and [V] we apply the formalism to the specific case of the Henon-Heiles 
hamiltonian. The numerical results exhibited in section [V| fully confirmed the expectations 
about the MNF as a tool for the study and characterization of chaotic motions. Moreover, 
they also point to some kind of extension of the region of convergence initially assumed for 
Moser's theorem. In fact, this issue is just being considered by the authors presently. 



Finally, in section [VI| we summarize the results and possible extensions of the present 
work. 
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II. MOSER'S NORMAL FORM 



It is essential that the hamiltonian be in the complexified form, for the method's imple- 
mentation: 



h(x 1 ,x 3 ,x 2 ,x A ) = XiXiX 3 + X 2 x 2 x A + H(£) x- 

e=3 



(1) 



Here, the origin is assumed to be a saddle point and we use the notation: £ = (£ 1 , £ 3 , £ 2 , £4), 
x = (xi, £3, x 2 , x 4 ), £ = £1 + £3 + £ 2 + £4 and x- = x^ 1 x e 3 x 2 2 x e 4 4 . The positions are x± and 
x 2 and the conjugate momenta are, respectively, x 3 and X4. The eigenvalues of the system's 
linear part are Ai = iuo and A2 = —A, with uo and A real (the other two being of course — iu 
and A). H(t) is the coefficient of x- and £ is its order. 

The usual noncomplexified form of the quadratic part h 2 of (|l|), around the saddle point, 



is 



^2(91, 92, Pi, P2) = ^(pl ~ ql) + ^(p? + ql) 



both sets of coordinates being linked by the symplectic transformation 



(2) 



(3) 



qi = 73(^1 + ^3) 
pi = 75(^1 + ^3) 
g 2 = -^{x 2 + x A ) 

P2 = 7= (~X 2 +X 4 ) . 

The Hamilton equations coming from ([!]) are: 

x x = Aixi + E^s^^^^^-Va 2 ^ 4 
x 3 = -X^-EZshHiDx^xfxi^x^ 
x 2 = \ 2 x 2 + Y.T= 3 hH{£)x^x e ix £ 2 2 x l r l 
x 4 = -A 2 x 4 -E^ 3 4^(I)4 1 4 3 4 2 "M 4 • 
Moser's theorem assures that there exists a near-identity polynomial coordinate trans 
formation, convergent in a neighbourhood of the origin, 



(4) 
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Xi = yi + J2 x ( z ^)y- 

1=2 



i = 1,2,3,4 , 

such that under it the original system (f|) takes the normalized form 



(5) 



(6) 



yi = F (vm , 2/22/4) yi 
m = -F{vm , 2/22/4) 2/3 

2/2 = G(yty 3 , 2/22/4) 2/2 
y A = -G(yxy 3 , 2/22/4) 2/4 , 

the functions F and G depending as indicated only on the products 2/12/3 and 2/22/4- 111 fact, 
the normal system (0) is also hamiltonian, with Hamilton function k(y) given by 



^(2/1,2/3,2/2,2/4) = A12/12/3 + A22/22/4 + 51 ^( 2 22l) (2/12/3) mi (2/22/4) 



1)12 



(7) 



m=2 



where K(2m) = i^(2m 1 ,2m 2 ) are the expansion's coefficients. 

It is straightforward to see that the two constants of motion of (|J) are just the products 



2/12/3 = ci 



(8) 



2/22/4 



c 2 



in terms of which the system can immediately be integrated: 



(9) 



(10) 



Vl {t) = yi (0) exp(tF) 
2/sW = 2/3(0) exp(-tF) 

y 2 ( t ) = y 2 (Q) exp(tG) 

y 4 (*) = 2/4(0) exp (-tG) , 

with F and G given by 

F{ yi y 3 , 2/22/4) = A 1 + E^=2^i^(2m)(2/i2/3) mi ~ 1 (2/22/4) m2 
G(2/i2/3, 2/22/4) = A 2 + E^=2^2^(2m)(2/i2/3) mi (2/22/4) m2 - 1 ■ 
It is important to note that the coordinate transformation |5] need not to be canonical, 
even though the transformed system is also hamiltonian. The cylindrical topology of the 



(11) 



flow can be made explicit. As the original coordinates qi,q2,Pi,P2 are real, we see from @ 
that 22 and 24 are real and x\ and 23 are complex, there existing the following relation 
between x\ and 23: 



-% X\ 



(12) 



The overline symbol above indicates the complex conjugate operation. Now, the first order 
truncation of ([5]) permits us to extend the above conclusions to the corresponding coordinates 
y, namely y 2 and 2/4 are real and y\ and y 3 are complex satisfying the relations 



2/3 



1 2/1 1 



m\ = p 



(13) 



2/12/3 



-%p 



(14) 



p being a new constant of motion. Then, it follows from ([H]) that F is pure imaginary and 
G is real. If we define (remembering that Ai = iuo and A2 = —A) 



Q = F/i = uj + 5u 
A = G = X + 5X , 



and 



5 a; 

sx 



-1 E~ =2 ^i^(2m) 
E^ =2 m 2J RT(2 m )(-zp 



2\m 1 (,\m 2 -l 



we can rewrite (|10D as 



(15) 



(16) 



(17) 



y x {t) = p exp[i(9 + nt)] 

y 3 (t) = pexp[-i{6 Q + * + m)] 

2/2 (*) = 2/20 exp [At] 

Viit) = y m exp [-At] , y 2 y\ = e , 

so that e is the second constant of motion. The convergence of the MNF guarantees that 
p and e are proper constants of the motion. The solution (|I7]) has an obvious cylindrical 
character: each given orbit, characterized by the phase #o> slides on the cylinder that is the 
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direct product of a circle with radius p in the plane 7/1, 7/3, with a hyperbola with parameter 
e in the plane y 2 , Va- Obviously, this topology is preserved in the original coordinates, due 
to the convergence of ([5]). For a linear flow generated just by h 2 in (0) all frequencies Q 
and velocities A degenerate respectively into uj and A, no matter what particular cylinder is 
considered. 

The suitable form of the solution (|17|) immediately reveals some important structures 
near the saddle point. Let us consider p = 0, in which case the motion is confined to the 
plane 7/2,7/4. For e = 0, we identify the coordinate axes 7/2 and 7/4 as just the stable and 
unstable orbits originating from the saddle point. For e 7^ 0, we have hyperbolae near the 
axes. 

Consider now the case p 7^ (see figure |l|). If y 2 = 7/4 = (e = 0), we have the family 
of unstable (circular) periodic orbits confined to the plane 7/1,7/3. Let r be one such orbit. 
If now only 7/4 = 0, then the motion is composed by r in that plane and the axis 7/2, thus 
generating from r the two unstable cylinders (U\ and U2). On the other hand, if only y 2 = 0, 
the composition of r with the axis 2/4 generates the pair of stable semicylinders (Si and S 2 ). 
Now, let us intersect any of the four semicylinders above with a plane transverse to it. As a 
consequence of the Poincare-Cartan theorem [HJ], the symplectic area or action of the closed 
irreducible curve, formed at the intersection, will equal that of r. In fact, this is valid for 
any irreducible curve over the semicylinder. On the other hand, it is a lagrangian surface 
i.e. all reducible curves on it have null action. Finally, if we compose a circular motion in 
the plane 7/1,7/3 with a hyperbola in the plane 7/2,7/4 (e 7^ 0), we do not have semicylinders 
but smooth cylinders near the saddle point, distinct from those associated to the orbits r. 

Of course, the "rectified" structures described above in the y coordinates, will appear 
distorted — yet preserving their topology — in the original coordinates. Moreover, the 
distorsions are such that, far away from the saddle point, the cylinders execute a very 
complicated tangle among themselves if the regime is chaotic ||15|| . Because they are all 
confined to a compact energy surface in phase space, transverse crossings of such cylinders 
eventually occur far from the origin. The existence of those transverse crossings are the 



source of chaotic motion in the continuum flow. So, an orbit at a transverse intersection 
of an unstable semicylinder with a stable one (Ui with 5*2 or U% with Si in the discussion 
above) is just a homoclinic orbit which tends to the respective orbit r as t — > ±00. In 
the case of the families of smooth cylinders in the neighbourhood of (Ui, S2) or (f^Si), 
the selfcrossings eventually lead to closed orbits with increasing periods, which accumulate 
on the homoclinic orbits. All these continuum structures underly the homoclinic tangle 
observed in a (convenient, yet generally nonplanar) Poincare section of the flow. 

Because of its convergence, the MNF formalism allows us to start from precise initial 
conditions near the saddle point. Moreover, the preceding discussion shows the necessity 
of the MNF for the general definition of a given cylinder in phase space. This permits us 



to extend the (un) stable cylinders beyond the vicinity of the origin [|15[ and proceed to an 
accurate quantitative study of the structures presented above. 



III. RECURRENCE RELATIONS 



To consistently determine the coefficients X(i,£) of (Q) and K(2m) of (|7|), we insert (^) 
in (Hj) and compare the result with the system directly obtained from the hamiltonian (0). 
We see that this task, simple in principle, is laborious if we are interested in an algorithm to 
compute the coefficents X(i, t) and K(2m), up to an arbitrary order. Anyway, the recurrence 
relations we obtain are the following: 

A(l,n)X(l,n)+B(l,n)K(2nx,2n 2 ) = 
(n 3 + 1) H(n h n 3 + 1, n 2 , n A ) + 

n-l 

9(n - 3) £ (£ 3 + 1) H{£ h h + 1, 4, Q 4 - 

t=2 

INT (n 12) 

0(n-4) £ K(2m)W^ , (18) 

m=2 
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A(3,n)X(3,n) + B(3,n) K(2n 3 ,2n 2 ) = 
-{ni + 1) + 1, n 3 , n 2 , n 4 ) - 

n-1 

©(« - 3) E (^i + !) + !> 4», 4, 4) 4 " 

INT(n/2) 

9(n-4) E ^(2ffi)W£ , (19) 

m=2 

A(2,n)X(2,n) + B(2,n) AT(2ni, 2n 2 ) = 
(n 4 + 1) n 3 , n 2 , n 4 + 1) + 

n-1 

0(n - 3) E + 1) #(4, 4, 4, 4 + 1) 4 " 

£=2 

INT(n/2) 

9(n-4) E , (20) 

m=2 

A(4,ra)X(4,n) + 5(4,71) K(2 ni ,2n 4 ) = 
-(n 2 + 1) 5(ni, n 3 , n 2 + 1, n 4 ) - 

n-1 

e(n-3)E (^ + 1) 5(4, 4, 4 + 1, h) 4 - 

£=2 

INT(n/2) 

6(n-4) E ^(2m)^4-n • (21) 



m=2 

The definitions involved here are: 



INT(x) = integer part of x , 

9(x) = 1 i/ x > or i/ x < , 

n) = Ai(ni - n 3 - 1) + A 2 (n 2 - n 4 ) 

^4(3, n) = Ai(ni - n 3 + 1) + A 2 (n 2 - n 4 ) 

^4(2, n) = Ai(ni - n 3 ) + A 2 (n 2 - n 4 - 1) 

^4(4, n) = Ai(ni - n 3 ) + A 2 (n 2 - n 4 + 1) , 



(22) 
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5(1, n) 


1 77-3 714 


5(3, n) 




5(2, n) 


= n 2 5X~ 1 


5(4, n) 


= -^4« +1 



(23) 



5 l . = l if i = j or if i ^ j 



Wf^ = m\{n\ — riz) X(i, m — mi + l,n 3 — m\ + l,n 2 — m 2 ,n± — m 2 )x 

X©(ni — nil + 1) 0(^3 — m i + 1) 0(^-2 — m 2) 0(^-4 — m 2 ) + 

m 2 (n 2 - n 4 ) X(i, n x — mi, n 3 - mi, n 2 - m 2 + 1, n 4 - m 2 + 1) x 

xO(rii — mi) Q(n 3 — mi) Q(n 2 — m 2 + 1) 0(n 4 — m 2 + 1) . (24) 

Given n and m, the coefficients Wf^ depend only on the X(i,k) of order 
k = n — 2m + 2. As the minimum value of m is 2, it follows that k < n — 2 and hence k < n. 
As the maximum value of m is INT(n/2), it also follows that k > 2. 

The coefficients Z| arise through the relations 

£ = My)]* 1 Mj/)]' 3 [x 2 (j/)]^ [x 4 (g)]* 

oo oo 

= [ yi + £ i) j/f 1 [y 3 + J2 X(S, jj yi-f x 

i=2 i=2 
oo oo 

x [y 2 + J2 X(2, k) jfl* [y A + £ X(4, m) j^]* 

fc=2 m=2 
oo 

= E 4^ • (25) 

n=£+l 

Obtaining the coefficients is the heaviest computational task of the whole algorithm. 
The reason, simply stated, is that we have to obtain a series resulting from multiplying four 
terms, each of them being itself a series powered to an arbitrary integer! Given an index 
vector I, the coefficients Z| depend only on the X(i,k) with k < n. To obtain the Z| 
up to order n = N, the series in ( plf) needs to be truncated at combined minimal orders 
which depend on N. If iV changes (increases), these series' minimal truncations also change, 
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modifying the values of various Z| of order n < N previously obtained. Then, for each 
increase of N, we need to compute all Z| again, since we do not know (at least at the 
present stage) what coefficients, coming from the previous step, will not be changed by the 
actual calculation. Nevertheless, the process evidently converges as N —>■ oo. 

Now, let us consider a subset of coefficients characterized by having some null indexes 
in n = (ni, ri3, ri2, 714). It is easy to see that this subset depends only on the X(i,k) with 
k having the same null structure as n. Then, the more null indexes the vector n has, the 
easier we can compute the corresponding subset from the smaller number of coefficients 
of the respective subset X(i,k). Now, let us return to the cylindrical structures discussed in 
the end of the last section. We need only the coefficients X(i, m, n 3 , 0, 0) or X(i, 0, 0, n 2 , n 4 ) 
to respectively describe, in the original coordinates, the normalized planes y% = = or 
2/1 = 2/3 = 0. Unfortunately, the coefficients X(i, ni, n 3 , 0, 0) and X(i, 0, 0, n 2 , n 4 ) depend 
by themselves on the remaining ones, through their relations with the coefficients K(2m). 
So, we find that the decoupling occurs only if just one index is nonnull, as for e. g. the 
subsets X(i, 0, 0, n 2 , 0) or X(i, 0, 0, 0, n 4 ). In particular, the later coefficients are just those 
we need to compute the (un) stable manifolds originated at the saddle point. In that case, 
the corresponding coefficients Z| can be obtained up to a much higher order than the full 
set, as we will see in section |V[ 



Let us check the consistency of the recurrence relations (pf) to (0). Each first member 
is a linear combination of the unknowns X(i,n) and K(2m), and the second members de- 
pend on the known hamiltonian coefficients H(£), on the X(i,n) through Z^ and Wf^ and 
explicitly on the K(2m). Nevertheless, the appearance of X(i,n) and K(2m) in the second 
members occurs at lower orders than in the corresponding first ones, as expected. Addition- 
ally, they have the following property: if A(i,n) ^ then B(i,n) = and reciprocally. This 
is due to the eigenvalues Ai and A 2 being independent over the reals. So, those relations 
naturally split into two cases: firstly, we have A(i,n) ^ and B(i,n) = 0. This case allows 
us to obtain the corresponding X(i,n) in a direct way. 

The case A(i,n) = and B(i,n) ^ can only occurs for odd orders n. This case 
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has redundancies so, in order to obtain all the K(2m), it is only necessary to work with 
either the pair of eqs. (|T8| ) and ( p0|) or (|19|) and (|21|). The coefficients X(i,n) that remain 



undetermined in this case are the following: X(l, ki + 1, ki, k 2 , k 2 ), X(3, k±, k x + 1, k 2 , k 2 ), 
X(2, ki, ki, k 2 + l, k 2 ) and X(4, k±, ki, k 2 , k 2 + l), with k± + k 2 = 1, 2, 3, ■ ■ •. In fact, the Moser 
theorem imposes some additional normalizing conditions to these "weak" coefficients, which 
can be simply met by requiring that all of them vanish. This will be assumed here. It should 
be noted that requiring the coordinate transformation to be canonical reduces this freedom 
and may not be compatible with the conditions of Moser's theorem. 

IV. PREPARATION OF THE HENON-HEILES HAMILTON! AN 

We apply the preceding formalism to the Henon-Heiles hamiltonian: 

h(Qi,Q2,Pi,P2) = ~(pl +Pl) + U(qi,q 2 ) , (26) 

U(qi,q 2 ) = \{ql + ql + H<l2-\ql) • (27) 

Here, qi and q 2 are the positions and p\ and p 2 are their respective momenta. 

Even though the motivation for the hamiltonian ( PB"D was astronomical, it soon became 
relevant on its own, due to both its simplicity and its dynamical richness [U,H,|12,|17| . For, 



if the energy h = E is suficiently small, the chaotic regions of the flow are vanishingly 
small. On the other hand, if the energy increases, the chaotic regions arbitrarily increase 
too. Additionally, (|26|) comes from the truncation of the Toda hamiltonian (see e. g. [|12|j), 
the later one being totally integrable for all energies. 

In figure we exhibit the equipotential curves of (|27D. The origin Po(<?i = 0, pi = 0, q 2 = 
0,p 2 = 0) is an eliptic (stable) equilibrium point with the eigenvalues ±i, while P\(0, 0, 1, 0), 
P 2 (-V / 3/2,0, -1/2,0) and P 3 ( v / 3/2,0, -1/2,0) are saddle (unstable) points, all having the 
same set of eigenvalues ±«v3, ±1. These saddle points are just as required by our method. 
The evident trigonal symmetry of figure |2| makes them mutually equivalent from a dynamical 
point of view. Then, there is just one MNF valid around P 1; P 2 and P3. 
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To put the hamiltonian (p6|) in the suitable form (|l|), we need to perform the following 
coordinate transformations for Pj , j = 1, 2, 3: 



pi 



( ~ \ 

Qi 
Pi 

\P2j 



+ P 



(28) 







( ~ \ 


Qi 




Z2 


Pi 


= 4- 


%1 




h 


Z A 


\P2) 




[z 3 ) 



(29) 



Zl 



Z3 
Z 2 

\~ Za J 

( \ 

Zl 

zz 
z 2 

\ Z4 J 



^ ^ 
^3 
10 
1 



1 

72 



J 

1 i 
i 1 
11 
0-11 



/ \ 

Zl 
Z3 
Z 2 

\ Za J 



\ ( \ 

X\ 

^3 
\ X4 J 



(30) 



(31) 



Here, Pj is the transpose of the respective saddle points given in the text and the matrices 
Aj are given by 



Ai 



^0 10 0^ 



10 
1 
10 



(32) 
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if the point is P±, and 



4 



cosOj sinOj 



\ 



cosOj sinOj 



\ 



J = 2,3 



(33) 





-sinOj cosOj 


—sinOj cosOj 

if the point is P 2 or P3, with 6 2 = ^ and 9 3 — |. All transformations made here are 
symplectic. 

In the intermediate z coordinates, the hamiltonian form is 



h(z) 



1 



x z 2 A - zl) + —{z\ + z 2 :) + —ziz 2 zl + - 



(34) 



The dissociation energy E = 1/6 arises explicitly here. The final hamiltonian form, in the 
x coordinates, is 

h(x 1 ,x 3 ,x 2 ,x i ) = i\Z3xix 3 - x 2 x 4: + ^(x 1 + ix 3 ) 2 (x2 + X4) -^(x 2 + x i ) 3 , (35) 

where we have dropped the constant energy term. We obtain the coefficients Ai, A2 and H(£) 
by identification of (p5|) with (JJ). We see that, for each Pj, we arrive at the same prepared 
hamiltonian form (|35| ) and hence we have a single MNF valid for all three Pj, as expected. 

V. RESULTS 



To present the results for the Henon-Heiles hamiltonian, we will use the z coordinates 
of the last section. All numerical propagations were made via an optimized fourth order 
Runge-Kutta code |18| and, for the Poincare sections, we use a trick by Henon [19|] for exact 
crossings of the surface sections. 

First of all, we compute the coefficients of the special form X(i, 0, 0, n 2 , 0) and 
X(i, 0, 0, 0, 714). As we saw in section |T|, this computation is much easier than for the 
full set X(i,n). So, with the same machinery and twice the real time for the general case, 
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we arrive only at the order 16 for the full set, against the order 50 for both subsets. The 
maximum order obtained for the coeficients K(2m) was also 16. 

Figures [5] and |] show both the exact (un) stable manifolds that originate at the saddle 
point and their approximations via the coordinate transformations z(y) (see fl31|) and (|5])), 
truncated at various orders up to the maximum of 50. We see that the exact curves meet 
smoothly in the bounded region. In fact, they pertain to a family of orbits entirely confined 
to the plane Z\ = = 0, as is easily seen from the form fl34f ) of the Henon-Heiles hamiltonian. 
To obtain the exact curves, we select one starting point near the origin at each axis t/2 or y^ 
(see ([H])) and then propagate them through the exact flow generated by h(z). This method 
will be referred to in the following as z-propagation, being always precise, provided we start 
from y near the origin. The z-propagation is semianalytical in the sense that we numerically 
evolve precise initial conditions that are analytically given. We note that for localizing the 
"correct" initial directions i.e. those given by the axes y2 and 2/4, the form of ( |P7D is crucial. 
On the other hand, the approximate curves in figures [| and ^ are just the axes yi and 
?/4, merely rewritten in the coordinates z via the series truncated at different orders. This 
last method will be referred to in the following as yz-translation and its precision is only 
guaranted by Moser's theorem if y remains in the neighbourhood of the origin. The yz- 
translation is analytical in the sense that it uses only algebraic calculations and properties 
of the normal form. 

The remarkable feature of both figures is that the analytic yz-translation clearly points 
to the convergence of the series (|5|) even far from the small region where Moser initially 
guaranteed it. This (and additional evidence in the following) strongly suggests the possi- 
bility of extending Moser's original convergence region. This issue is being considered in our 
present research. We also note that the region of accurate approximation ceases more and 
more abruptly as the truncation order become higher. 

The next task is to obtain the unstable periodic orbits near the saddle point. They 
are entirely confined to the plane yi = y^ = (e = 0) of the normalized system, the fam- 
ily parameter being p (see (|I"7|)). We need in principle only the coefficients of the form 
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X(i, ni, ris, 0, 0) to compute them. However, we saw in section [H| that those and all remain- 
ing coefficients must jointly be found. We obtained the full set X(i, n) in this work until 
the order 16 and from now on this will be ever the case. 

Figures |5| and [| show the projections of three such unstable closed orbits in the planes 
z 2 = z^ = and z\ = z% = 0, respectively. The values of p caracterizing each orbit are 
shown in the figures. For sufficiently small values of p we are always in the convergence 
region for 2/2 = 2/4 = and both methods i.e. the z-propagation and yz-translation, generate 
closed curves which are indistinguishable in practice. Any orbit of the family closes once in 
figure |5| for every two returns in figure || We also note that the motion is microscopic or 
"residual" in the second figure (in the plane z\ = Z3 = 0) relatively to the first one (in the 
plane z 2 = z^ = 0). These features are easily understood if we consider y 2 = 2/4 = in ([|) 



and (|31|) and keep only contributions from the £ = 2 terms in (g). 

We saw in section [HI] that there exist a cylindrical flow associated to any of the unstable 
periodic orbits above. Four semicylinders emanate from each of them. Let r be one such 
orbit and consider, for example, only the pair of (un)stable semicylinders which directly 
evolve from r to the region z 2 < 0. This evolution occurs in an intrincate manner in the 
energy shell and, because the dynamics are chaotic, the semicylinders mutually intersect 
transversally. At each of these intersections there exists a homoclinic orbit (HO) which 
tends to r as t — > ±00. 

We consider in this work only the first of the intersections above, in which we have the 
primary HOs. A suitable Poincare surface to observe that is one that is transverse to the 
cylinders themselves, i.e., one that intersects each cylinder in an irreducible (closed) curve. 
The symmetry of figures |3| and ^ and the continuity of the cylinders with respect to the 
parameter p suggest that one such surface is the plane z^ = 0. On the other hand, figure |6| 
shows that this plane also intersects the corresponding orbit r ( and hence the semicylinders 
themselves) in the region z 2 > 0. In fact, it is seen both qualitatively and numerically that 
in the region z 2 > the cylinders do not mutually intersect at all, except obviously at r 
itself. Moreover, the plane z± = generates only reducible (open) curves on the cylinders in 
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that region i.e. it is longitudinal to the cylinders there. For sufficiently small p the orbit r 
is near the origin, so we avoid numerical troubles, without losing precision, by starting on 
the cylinders at some small z 2 < 0. 

We obtain the Poincare section z^ = in the region Z2 < by three methods. The first 
one is the z-propagation of the following initial conditions: a sufficiently small (constant) 
height ?/4 = 7/40 determines a straight section on the semicylinder characterized by 7/2 = 
(hence e = 0) and a given p (see (|17|)). Now, all orbits through the section are given by 

?/i = p exp \i9] 

[ J (36) 

7/3 = -ip exp [—i0] , 9 G [0, 2tt) , 
where the phase 9 identifies each orbit on it. Obviously, we translate all initial points above 
from 7/ to z, before propagating them by the exact flow. The propagation of the other 
semicylinder, characterized by 7/4 = and the same p, is made in a similar way, except 
that the straight section is now determined by the small height 7/2 = 7/20- In all cases 
presented here it proved to be sufficient to use 7/20 = V40 — 1.0 x 10~ 2 . We z-propagated 
both semicylinders until the condition z 4 = was reached by the first time in the region 
z 2 < 0. 

Figure [1] shows a typical Poincare section as described above, namely that for p = 
1.0 x 10 -3 . Four primary HOs are exhibited there. The highly symmetric features of this 
figure (and all succeeding ones) must be viewed as reflecting the intrinsic symmetries of the 
Henon-Heiles potential. So, they do not depend on the specific value of p, which works 
somewhat as a scaling factor in all them. In particular, we see that the HOs are always on 
the indicated in figure [7|. Each of these HOs emanates from the corresponding orbit 

t along the unstable cylinder and returns to r along the stable one. The remaining orbits 



do not exchange cylinders, at least at the first crossing. Figures p|, |9|, |10| and [TT| show, for 
p = 1.0 x 10~ 3 , the evolution of a typical orbit on each cylinder until the Poincare section 
was reached. On the other hand, figures 0, [13] and [L4] show a primary HO associated to the 
orbit r with p = 2.6 x 10 -3 , just that which corresponds to the HO localized on the left (i.e. 



at Z3 < 0) in figure 0. Figure |H| is an amplification of figure 0, revealing the microstructure 
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of the HO near r. We note that two typical orbits on distinct semicylinders do not have the 
same coordinates when they reach, for the first time, the condition z^ = 0, contrary to what 
occurs with the two branches of a primary homoclinic curve. 

We also consider two analytical methods for the determination of the same Poincare 
sections above. One of them is the yz-translation already presented. For this purpose we 
have merely to rewrite the condition z^ = in the y coordinates, 

N 

yf +1 + ^[X(2,n 1 ,n3,0,n 4 )-X(4,n 1 ,n3,0,n4M 1 l/3 3 l/f +n4 = , (37) 

n=2 

for the semicylinder for which y 2 — and 

N 

^ +1 + ^[X(2,n 1 ,n 3 ,n 2 ,0)-X(4,n 1 ,n 3 ,n 2 ,0)K 1 l/3 3 l/2 A ' +n2 = , (38) 

n=2 

for the one with y± = 0. Here, N = 16 and y\ and y% are fixed through (|36"D for each value of 9. 
The equations (|37| ) and fl38|) thus determine for each phase/orbit 9 what is the corresponding 
height 2/4 or y 2 where it crosses the plane = 0. As they are polynomials of degree 2N in 
?/4 and y 2 respectively, we have to select in each case, among the 2N complex roots, the only 
one that is real and verifies z 2 < 0. To find the roots, we use codes described in |T8T], based 



on the Laguerre method. Finally, we come back to the results in the z_ coordinates. 

The second analytical method is essentially equivalent to the yz-translation and consists 
of evolving in time the solutions (0) themselves until the condition z 4 = be verified. We 
will refer to it as the y-propagation. Even though the y-propagation be trivial, the com- 
parison of both analytical methods is an important test for the numerical consistency of 
both sets of coefficients X(i,n) and K(2m) that we obtained, as the yz-translation involves 
only the coefficients X(i,n). In fact, the figures for the Poincare sections obtained through 
both analytical methods are indistinguishable in practice, which confirms the mutual consis- 
tency of those two sets of coefficients. On the other hand, in figure [15| we superimpose (for 
p = 1.0 x 1CT 3 ) the analytical and the semianalytical results, the latter already displayed in 
figure 0. We conclude that the analytical method is able to describe all qualitative aspects, 
such as symmetries, and the existence of the transverse cylindrical intersections and HOs 



- and hence the chaotic nature of the motion. Moreover, we see that even its numerical 
efficiency is surprising if we realize that we are far from the region of guaranteed conver- 
gence. Considering that here the MNF has only been expanded to order 16, we are faced 
again with an indication that Moser's theorem can be extended. 

Now we are ready to search for the closed orbits with long periods, which accumulate 
on the HOs. For this purpose, we have to consider the smooth selfintersecting cylinders for 
which not only p ^ but also e ^ (see section |TTT|) . Due to the continuity of the cylindrical 



structure with respect to the parameters p and e, we know that the transverse cylindrical 
selfcrossings must also occur at the plane z^ = in the region z% < 0. Then, all the previous 
discussion about both the semianalytical and analytical methods are still in order here, the 
only changes being now that the relation y 2 Z/4 — e eliminates both conditions (|37|) and (p8|) 
in favour of just one, e. g. 

N 

^ +1 - e ^- 1 + ^[X(2,n)-X(4,n)]< 1 2/3 3 ^ +n2 - n4 = , (39) 

n=2 



with yi and y 3 still given by 

Figure |H] shows the Poincare section above, obtained through the semianalytical method 
for the cylinder with parameters p = e = 1.0 x 1CT 3 . It exhibits four orbits in the first 



autointersection of the cylinder, with the same symmetries already seen. Figures [T7| and [L8 
show a typical orbit on that cylinder, from the time it crosses the plane z^ = through the 
branch of the cylinder which brings the orbit near to the origin (the S-branch), until the time 
it crosses that plane again, coming through the branch which takes it away from the origin 
(U-branch). We compared with each other the two analytical methods (yz-translation and 
y-propagation), confirming also in this case the underlying numerical consistency of both 
sets of coefficients X(i,n) and if (2m). In figure [19] we superimpose the analytical results 
with those displayed in figure [L6|. 

What happens to the orbits in the selfintersections of figure |16| ? We easily discover by 
numerical computation that, they do not close in general for any specific values of p and 
e — at least at the first intersection. In fact, as the S and U-branches of the cylinder 
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met smoothly near the origin, a more complex evolution can occur to any orbit before it 
eventually closes after a certain number of cylindrical selfmtersections. We consider here 
the accumulation of periodic orbits only on the primary HOs, so we will limit the search of 
closed orbits just at the first crossing. For this reason, we fix p at some value and vary e, thus 
verifying for each e whether the four primary autointersections close on themselves or not. 
They in fact do so for a set of well defined discrete values of e. In particular, the four curves 
close simultaneously on themselves when they do so, generating four distinct periodic orbits 
of the same period. We display in table | the numerical values of energy and period for the 
case p = 2.6 x 10~ 3 . We present there only those orbits at the crossing localized on the left 
(with z 3 < 0) in figure [16|, while the corresponding primary HO of interest lies at the same 
crossing, as in figure [7[ P01,...,P04 are the periodic orbits we found, E is the energy and T 
is the period of each orbit. For the HO the finite numerical value for T is also shown. Due 
to the high numerical instability associated to the chaotic regime, the orbit are quickly lost 
after that time. We note the fast convergence of the values of E and also that the periods 
of the POs are (nearly) multiple of that of r. On the other hand, the accumulation process 
of P01,...,P04 onto the corresponding HO, which in turn tends to r as t — > ±00, is clearly 
shown in figures 120, BT] and 123. These figures are typical of the behaviour of the periodic 



orbits for all other values of p. 

Finally, an independent test was made to assure that the periodic orbits above are 
computed with precision by the present method, i.e. we ask whether the truncation at 
N = 16 of the series (||) is sufficient or not to yield precise results. For this sake we use 
a code |20| already known from the literature for computing periodic orbits in hamiltonian 



flows. By starting from a guessed periodic orbit, the code searches for an exactly closed orbit 
in its vicinity (in terms of either the energy or the period). That code is based on rewriting 
the linear finite difference integration of the flow in terms of the monodromy matrix and its 
main feature is the very fast convergence rate. 

So, we enter as guessed periodic orbits, the closed ones we have obtained above and 
compare each of them with the corresponding exact curve we get through the code. It turns 
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out that both orbits are indistinguishable in all cases we considered. For example, the exact 



curves fit completely those shown in figures |2(] to |22|, at those scales. Hence, we conclude 
that the truncation at iV = 16 we have used is sufficient to support the semianalytical results 
presented here. 



VI. CONCLUSIONS 

We have applied the theory of Moser's convergent normal forms to a two-degree of 
freedom hamiltonian flow near a saddle point. It turns out that the cylindrical topology of 
the flow is promptly revealed in the normal coordinates. The chaotic behaviour was studied 
directly in the continuous flow underlying the homoclinic tangle in a Poincare section. We 
applied the method to the Henon-Heiles hamiltonian. This allowed to compute structures 
related to chaotic motion in the full phase space rather than just a Poincaree section: the 
(un)stable manilfolds originated at the saddle point, the unstable periodic orbits with the 
cylinders they generate and the homoclinic orbits associated to them, and also the closed 
orbits with long periods which accumulate on the homoclinic ones. 

We developed two methods for obtaining those structures. The first one is semianalytical 
in the sense that we numerically propagate initial conditions that we were able to find only 
through the normal form. The second method (in two versions — which are important to 
test the numerical consistency of all coefficients obtained) is entirely analytical, in the sense 
that it uses only algebraic calculations and properties of the normal form. The calculations 
are truncated at the order 16 and compared with a known numerical method for searching 
exact periodic orbits in a two-degree hamiltonian flow. Our semianalytical results proved 
accurate, thus showing the sufficiency of the truncation order that we have used. On the 
other hand, the analytical method was surprisingly good even far away from the saddle point 
and it was certainly able to exhibit all the qualitative aspects of the chaotic motion there. 

It was also possible to compute analytically, as a particular case, the (un) stable manifolds 
originated at the saddle point up to the order 50, obtaining much better numerical agreement 
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for the analytical method. Both facts indicate the convergence of the normal form beyond 
the "microscopic" region originally established by the Moser's theorem. We are presently 
working out this extension. 

Despite the mathematical interest in ascertaining the convergence of Moser's normal form 
as far as the first intersection of the stable and unstable manifolds, there is no doubt that, 
even for a simple polynomial potential as the Henon-Heiles, the convergence will be very 
slow. In view of the huge difficulty in increasing the order of the normal form, it seems that 
the semianalytical method that we have used to obtain chaotic structures is to be preferred. 
This relies only on a relatively small truncation of the normal form near the saddle point, 
which produces accurate cylinders that may be extended by numerical integration. 
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TABLES 





e 


E 


T 


P01 


5.836022 x 1(T 5 


0.1666200153 


14.026 


P02 


4.1216 x 1(T 8 


0.166678334123 


21.381 


P03 


2.9118 x KT 11 


0.166678375310 


28.536 


P04 


2.057124 x 10~ 14 


0.166678375339 


35.791 


HO 





0.166678375339 


(oo) 83.795 


r 





0.166678375339 


3.628 



TABLE | 
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TABLE CAPTIONS 
TABLE I. Parameters of some orbits discussed in the text with p = 2.6 x 1CT 3 . 
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FIG. 1. Projections in the planes yi 
near a saddle point. 



FIGURES 

= y4 = (a) and y\ = 7/3 = (b) of a hamiltonian flow 



FIG. 2. Equipotential curves of the Henon-Heiles potential shown in eq. (p7|). 

FIG. 3. Exact (un)stable manifolds that originate at the saddle point and their MNF approxi- 
mations at various even truncation orders 2, • • -, 50. The exact curves meet smoothly at Z4 = —1.5 
and a particle in the loop runs in the clockwise direction. Note the evidence of convergence of the 
MNF far from the origin. 

FIG. 4. The same as in the figure || except that the truncations are made at the various odd 
orders shown. Note in this case the "out" direction of the divergence of the approximate curves 
from the exact one. 

FIG. 5. Projection in the plane z^ = z^ = of three unstable periodic orbits (for which 
2/2 = 2/4 = and hence e = 0) near the saddle point, for the values of p indicated in the figure. 

FIG. 6. Projection in the plane z\ = 2:3 = of the orbits shown in the figure ||. Note that two 
loops here corresponds to one loop in the other projection. Note also the microscopic motion here, 
relatively to that in the plane Z2 = z^ = 0. 

FIG. 7. Poincare section at Z4 = (in the region Z2 < 0) of the (un)stable semicylinders (for 
which 2/22/4 = e = 0), originated at the orbit r with p = 1.0 x 10 -3 . Note the four primary 
homoclinic orbits at the intersections on the axes. 

FIG. 8. Projection of the (time-reversed) evolution of a typical orbit on the stable semicylinder, 
from a small initial condition already in the region Z2 < 0, until the condition Z4 = was reached. 
The parameters are p = 1.0 X 10~ 3 and y22/4 = e = 0. 



FIG. 9. The other projection of the orbit shown in the figure 
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FIG. 10. Evolution of a typical orbit on the unstable semicylinder in the same conditions of 
the figure ||. 



FIG. 11. The other projection of the case shown in the figure 10 



FIG. 12. Projection of a primary homoclinic orbit (the HO indicated in the text), associated 
to the orbit r with p = 2.6 x 10~ 3 . The computed time was T = 83.795, much longer than the 
period T = 3.628 of t, and its energy E = 0.166678375339 is the same of r, as expected (see table 



FIG. 13. Amplification of the figure 12 showing the microstructure of the HO near the orbit r. 



FIG. 14. The other projection of the homoclinic orbit shown in figure ff2l compared with the 
corresponding orbit r. 

FIG. 15. Superimposition of the precise semianalytical results, already displayed in the figure [?], 
to that obtained by the analytical method. Note the qualitative agreement and even the surprising 
numerical proximity between them, although we are here very distant of the origin. This is an 
additional evidence for the convergence of the Moser normal form beyond that neighbourhood. 

FIG. 16. Precise semianalytical obtaining of the Poincare section at = (in the region 
Z2 < 0) of the (smooth) cylinder for which p = e = 1.0 x 10~ 3 . The resemblance with the figure fj] 
is due to the intrinsic simmetries of the Henon-Heiles potential. 

FIG. 17. Projection of a typical orbit on the (smooth) cylinder for which p = e = 1.0 x 10 -3 , 
from a small initial condition already in the region Z2 < 0, until the condition z<± = was reached 
in both time directions. 



FIG. 18. The other projection of the orbit shown in the figure |17[ Note that it closes on itself 
on the plane shown in that figure but not here. 



28 



FIG. 19. Superimposition of the results already displayed in the figure 16, to that obtained by 
the analytical method. The comments made in figure about the comparison of the semianalyt- 
ical and analytical methods, are in order here. 

FIG. 20. Superimposition of the orbits of the table |, showing the accumulation of the periodic 
orbits (POl to P04 in the table) onto the corresponding homoclinic orbit (HO), associated to the 
r orbit with p = 2.6 x 1CP 3 . There is no difference among them far from the origin at this scale. 



FIG. 21. Amplification of the figure |20| near the origin. The accumulation process is much more 
slow here, on the contrary that occurs in the far distant region. The periodic orbits POl and P02 
of the table | does not reach this region. 

FIG. 22. The same accumulation process of figure ^ projected now onto the other plane. These 
results are typical of all others values of p we have considered. 
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